GLOBAL WELL-POSEDNESS OF THE THREE-DIMENSIONAL VISCOUS AND 
INVISCID SIMPLIFIED BARDINA TURBULENCE MODELS 
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Abstract. In this paper we present analytical studies of three-dimensional viscous and inviscid simpli- 
fied Bardina turbulence models with periodic boundary conditions. The global existence and uniqueness 
of weak solutions to the viscous model has already been established by Layton and Lcwandowski. How- 
ever, we prove here the global well-posedness of this model for weaker initial conditions. We also establish 
an upper bound to the dimension of its global attractor and identify this dimension with the number 
of degrees of freedom for this model. We show that the number of degrees of freedom of the long-time 
dynamics of the solution is of the order of (L/l d ) 12 / 5 , where L is the size of the periodic box and l d is 
the dissipation length scale- believed and defined to be the smallest length scale actively participating 
in the dynamics of the flow. This upper bound estimate is smaller than those established for Navier- 
Stokes-a, Clark-a and Modified-Leray-a turbulence models which are of the order (L/i^) 3 . Finally, we 
establish the global existence and uniqueness of weak solutions to the inviscid model. This result has 
an important application in computational fluid dynamics when the inviscid simplified Bardina model 
is considered as a regularizing model of the three-dimensional Euler equations. 
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1. Introduction 

Let us denote by v(x,t) = (vi(x,t),V2(x,t),Vs(x,t)) the velocity field of an incompressible fluid and 
p(x, t) its pressure. The three-dimensional (3D) Navier-Stokes equations (NSE) 

d t v — vAv + V • (v ® v) = — Vp + /, 

V-u = 0, (1) 

v(x,0) = v m (x), 

governs the dynamics of homogeneous incompressible fluid flows, where f(x) = (fx(x), f2(x), fa(x)) is 
the body force assumed, for simplicity to be time independent. The existing mathematical theory and 
techniques are not yet sufficient to prove the global well-posedness of the 3D NSE. Researchers who are 
investigating this question have incorporated the use of computers to analyze the dynamics of turbu- 
lent flows by studying the direct numerical simulation (DNS) of these equations. However, this is still 
a prohibitively expensive task to perform even with the most technologically advanced state-of-the-art 
computing resources. Tracking the pointwise flow values by numerical simulation for large Reynolds 
number is not only difficult but also, in some cases, disputable due to sensitivity of numerical solutions 
to perturbation errors in the data and the limitations of reliable numerical resolution. In many practical 
applications, knowing the mean characteristics of the flow by averaging techniques is sufficient. However, 
averaging the nonlinear term in NSE leads to the well-known closure problem. 
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To be more precise, if v denotes the filtered/averaged velocity field then the Reynolds averaged NSE 
(RANS) 

d t v - vAv + V • (vWv) = -Vp + f, 

V-w = 0, ^ 

where 



V • (v ® v) = V • (v ® v) + V • 72-(v, u), 
7c(i>, u) = v ® t> — w ® « 

is not closed. The quantity lZ(v,v) is known as the Reynolds stress tensor. The RANS system of equa- 
tions contains the unknown quantity v = v — v, which represents the fluctuation around the filtered 
velocity v. The equation in (23) is not closed because we cannot write it in terms of v alone. The main 
essence of turbulence modeling is to derive simplified, reliable and computationally realizable closure 
models. 

In 1980, Bardina et al. [3] suggested a particular closure model by approximating the Reynolds stress 
tensor by 

lZ(v, v) « v ® v — v v. (4) 
In [31], Layton and Lewandowski considered a simpler approximation of the Reynolds stress tensor, given 

by 

TZ(v, v) ps v <8> v — v Cg) v. (5) 
This is equivalent form to the approximation 

V • (vWv) w V • (v (git;). (6) 
Hence, Layton and Lewandowski studied the following sub-grid scale turbulence model: 



wt — vAw + V • (w CE3 w) = — Vq + f, 

V • w = 0, (7) 
w(x,0) = v (x), 

where they denoted (w,q), the approximation to (v,p). In this paper we will call this particular model 
the simplified Bardina model. Similar to the alpha models [6, 7, 8, 16, 10, 26], Layton and Lewandowski 
[31] used the smoothing kernel associated with the Helmholtz operator (J — a 2 A) -1 . That is, if v 
denotes the unfiltered velocity and u denotes the smoothed filtered velocity, then we have the relationship 
v = u — a 2 Au. For abstract mathematical study, one can define a more general smoothing kernels, which 
gives a different relationship between u and v (see, e.g., [4], [36]). In this paper, we will keep the 
same exact smoothing operator. There is a very important reason behind the choice of this particular 
smoothing kernel in our mathematical studies. The reason can be traced back from the early study 
of 3D Navier-Stokes-a (NS-a) turbulence model (also known as the viscous Camassa-Holm equations 
(VCHE) and Lagrangian averaged Navier-Stokes-a (LANS-a) model ). The explicit analytical steady 
state solutions to the NS-a model were found to compare successfully with empirical data for mean 
velocity and Reynolds stresses for turbulent flows in channels and pipes for wide range of Reynolds 
numbers (see, e.g., [6, 7, 8]). It was, in fact, this important finding, which led the authors of [6, 7, 8] to 
suggest that the 3D NS-a model be used as closure model for the Reynolds averaged equations (RANS). 
Under this particular relationship v = u — a 2 Au between u and v, the other alpha models reduced under 
the channel and pipe symmetry yield exactly the same equations, up to a modified pressure, to the system 
of equations for the NS-a model restricted to this symmetry. Hence, the explicit steady state solution 
to these equations will match the experimental data as well. This is one important property shared by 
all the alpha models. In particular, the simplified Bardina model enjoys this important property as well. 
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A more detailed discussion of this will be presented in section 3. With this at hand, we can rewrite the 
simplified Bardina model (7) as 

d t v — uAv + (u ■ V)u = — Vp + /, 
V • u = V • v = 0, 

v = u — a 2 Ait, (8) 
u(x,0) =u m (x) 
u and v are periodic, with periodic box Q = [0, 2ttL] 3 

Notice that consistent with all the other alpha models, the above system is the Navier-Stokes system of 
equations when a = 0, i.e. u — v. We have rewritten equation (7) in the particular form (8) in order 
to coordinate its similarity with the family of alpha models [6, 7, 8, 10, 11, 16, 23, 26]. In this form, 
when compared to the other alpha sub-grid scale turbulence models, the main difference, namely in the 
bilinear term, can be distinguished easily. 

Moreover, we note that, in addition to the remarkable match, in the channels and pipes, of explicit 
analytical steady state solutions of the alpha models to the experimental data the validity of the first 
alpha model, the NS-a model, as a subgrid scale turbulence model was also tested numerically in [9] and 
[35]. In the numerical simulation of the 3D NS-a model, the authors of [9], [19], [20] and [35] showed that 
the large scale (to be more specific, those scales of motion bigger than the length scale a) features of a 
turbulent flow is captured. Then, for scales of motion smaller than the length scale a, the energy spectra 
decays faster in comparison to that of NSE. This numerical observation has been justified analytically 
in [17]. In direct numerical simulation, the fast decay of the energy spectra for scales of motion smaller 
than the supplied filter length represents reduced grid requirements in simulating a flow. The numerical 
study of [9] gives the same results. The same results hold as well in the study of the Leray-a model in 
[10] and [19]. 

This paper is arranged as follows. In section 2 we fix some notations and define the functional setting. 
In section 3 we discuss in further details why we chose the particular smoothing kernel and justify the 
use of the simplified Bardina model as a closure model to the RANS. In section 4 we will re-establish the 
global existence and uniqueness of weak solutions of equation (8) subject to periodic boundary conditions. 
We will re-establish this result requiring a weaker initial condition than those required in [31]. In section 
5 we also provide an upper bound to the dimension of its global attractor. We then relate this upper 
bound to the number of degrees of freedom of the long-time dynamics of the solutions to this model. Our 
results show that the number of degrees of freedom for this model is proportional to (L/ld) 12 ^ 5 - This 
estimate is much smaller compared to those established for the 3D Clark-a [11], 3D NS-a model [16] 
and the 3D Modified-Leray-a model [26] which are of the order (L/ld) 3 - The smaller estimate on the 
number of degrees of freedom for the simplified Bardina model is expected since it has a milder nonlinear 
term than the 3D Clark-a model, 3D NS-a model and the 3D Modified-Leray-a model. Notice, however, 
that we have excluded the Leray-a model in our comparison above. For the Leray-a model, the estimate 
for its number of degrees of freedom is of the order (L/ld) 12 / 7 as shown in [10]. The power (12/7) is 
smaller than the power (12/5) of our estimate on the simplified Bardina model even though we have 
here a smoother nonlinear term u ■ Vu compared to nonlinear term u ■ Vv of the Leray-a. One reason 
for this is that the energy dissipation length scale Id for the Leray-a model is different from the Id of 
the simplified Bardina, 3D NS-a, Clark-a and Modified Leray-a model. For the Leray-a model, the 
dissipation length scale Id is based on the time average of the H 3 ~ norm of u. On the other hand, the Id 
of the simplified Bardina, 3D NS-a, Clark-a and Modified Leray-a model is based on the time average 
of the H 2 — norm of u. Recently, Holm and Gibbon [21] produced an interpretation of the dimension of 
the global attractor in terms of the Reynolds number. This global interpretation can assist in making 
across the broad comparison between the various alpha models. In particular, by following their work 
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one would be able to show that the dimension of the global attractor for the simplified Bardina model is 
much smaller than that of the NS-a model, but larger than that of the Leray-a. 

For completeness, in section 6 we also include in our study the energy spectra of the simplified Bardina 
model. Although the dimension of the global attractor for the simplified Bardina model is smaller in 
comparison to those established for the 3D NS-a model and the 3D Modified-Leray-a model, we found 
that the spectral slopes for the energy spectra for the simplified Bardina model is the same to that of 
3D Clark-a, 3D NS-a model and the 3D Modificd-Leray-a established in [11, 16, 26] respectively. 

In the last section we prove the global existence and uniqueness of the inviscid simplified Bardina 
model. This result has important consequences in computational fluid dynamics when the inviscid 
simplified Bardina model is considered as a regularizing model of the 3D Euler equations. This is 
because the inviscid simplified Bardina is globally well-posed model that approximates the 3D Euler 
equations without adding any hyperviscous regularizing terms. In particular, we propose the inviscid 
simplified Bardina model as a tool for testing claims about the formation of a finite time singularity in 
the 3D Euler equations (see, e.g., [25], [28] and references therein). 



2. Functional Setting and Preliminaries 

Let 0. = [0,2ttL] 3 . The simplified Bard ina turbulence model (8) of viscous incompressible flows, 
subject to periodic boundary condition, with basic domain f2, is written in expanded form: 

d t (u - a 2 Aw) - vA(u - a 2 Aw) + (u ■ V)it = -Vp + /, 

V •« = (), (9) 
u(x,0) = u in (x), 

where, u represents the unknown "filtered" fluid velocity vector, and p is the unknown "filtered" pressure 
scalar; v > is the constant kinematic viscosity, a > is a length scale parameter which represents the 
width of the filter. The function / is a given body forcing assumed, for the simplicity of our presentation, 
to be time independent and with mean zero, that is f(x)dx = 0, and u m is the given initial velocity 
also assumed to have zero mean and hence the solutions u and v as well. 

Next, we introduce some preliminary background material following the usual notation used in the 
context of the mathematical theory of Navier-Stokes equations (NSE) (see, e.g., [13, 44, 45]). 

(i) We denote by L p and H m the usual Lebesgue and Sobolev spaces, respectively. And we denote 
by | • | and (•, •) the L 2 — norm and L 2 — inner product, respectively. 

(ii) Let T be the set of all vector trigonometric polynomials with periodic domain Q. We then set 

V = jc£ G T : V • <j> = and^ dx = oj . 

We set H and V to be the closures of V in L 2 and H 1 , respectively. We also note that by Rcllich 
lemma (see, e.g., [1]) we have the V is compactly embedded in H. 

(iii) We denote by P a : L 2 — ► H the Hclmholtz-Leray orthogonal projection operator, and by 
A = — P a A the Stokes operator subject to periodic boundary condition with domain D{A) = 
(H 2 (il)) 3 n V. We note that in the space-periodic case, 

Au = -P a Au = -Au, for all u G D(A). 

The operator A" 1 is a self-adjoint positive definite compact operator from H into H. (cf. 

[13, 44]). We denote by < L~ 2 = Ai < A2 < the eigenvalues of A, repeated according 

to their multiplicities. It is well known that in three dimensions , the eigenvalues of the operator 
A satisfy the Weyl's type formula (see, e.g., [2, 13, 40, 45]) namely, there exists a dimensionless 
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constant co > such that 

J —<^<c j 2/3 , forj = l,2,...,. (10) 

We also observe that, D{A n / 2 ) = {H n {VL)f n V. 

(iv) We recall the following three-dimensional interpolation and Sobolev inequalities (see, e.g., [1] 
and [13]): 

Ml^cU^H^, and 
\\<t>h°<c\\cf>\\m, for every cpeH l (n). 
Also, recall the Agmon's inequality (see, e.g., [2, 13]): 

H0IU- ^cUW^UW 1 ^, for every (f> £ H 2 (Q). (12) 

Hereafter c will denote a generic dimensionless constant. 

(v) For wi,W2 G V, we define the bilinear form 

B(w 1 ,w 2 ) = P a ((w 1 -V)w 2 ). (13) 
In the following lemma, we will list certain relevant inequalities and properties of B (see [13, 44]). 

Lemma 1. The bilinear form B defined in (13) satisfies the following: 

(i) B can be extended as a continuous map B : V x V — > V , where V is the dual space of V . In 
particular, for every W\ , u>2 , w 3 S V, the bilinear form B satisfies the following inequalities: 

(B( Wl ,w 2 ),w 3 ) v , | < cKi^iKir^lMHKII, (14) 

(B(w 1 ,w 2 ),w 3 ) v , | < c|K|||K||| W3 r /2 |KH 1/2 . (15) 
Moreover, for every wi,w 2 ,ws € V , we have 

(B(wi,w 2 ),w 3 ) v , = - {B(w 1 ,w 3 ),w 2 ) v , . (16) 

And in particular, 

(B(w 1 ,w 2 ),w 2 ) v , = 0. (17) 

(ii) For wi s V and w 3 G D{A), we have 

{(B(wi,w 1 ),w 3 )v') | = | ((B(w uWl ),w 3 ) D{A y) I < A^V^alKIIKH, (18) 
where D(A)' is the dual space of D(A). 

Using the bilinear form B and the linear operator A, the sytems in (8) and (9) is equivalent to the 
functional differential equation 

dv 

— + vAv + B(u, u) = f, 

v = u + a 2 Au, ( 19 ) 

v(0)=v m = u m + a 2 Au m , 

Definition 2. (Weak Solution) Let f e H, u(0) = u m e V, and T > 0. A function u € C([0, T];V) n 
du 

L 2 ([0, T]; D(A)) with — G L 2 ([0,T]; H) is said to be a weak solution to (19) in the interval [0,T] if it 
dt 

satisfies the following: 

\^^ w ) +v(Av,w) D{Ay + (B(u,u),w) = (f,w), (20) 

\ at I D(A)> 
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for every w G D(A). Here, the equation (20) is understood in the following sense: 
For almost everywhere to,t € [0,T] we have 

(v(t),w) v ,- (v(t ),w) v ,+v j (v,Aw)+ f (B(u(s),u(s)),w)ds= f (f,w)ds. (21) 

3. The simplified Bardina model as a turbulence closure model 

As we mentioned earlier, one important characteristic shared by all the alpha models is the particular 
kernel used to give the relation between the smoothed velocity u and unsmoothed velocity v. This par- 
ticular choice of smoothing kernel gives the important result that under the pipe and channel symmetry, 
the reduced equation of all the other alpha models takes the form of the reduced of NS-a under the 
same symmetry, up to modified pressure. As a result, the explicit analytical steady state solutions to 
these equations will resemble the explicit analytical steady state solutions of the NS-a. In this way, the 
excellent match of explicit analytical steady state solutions of NS-a to experimental data in the channel 
and pipe symmetry for wide range of Reynolds number ([6, 7, 8]) is also inherited by these models. In 
this section, we consider the simplified Bardina model as a closure to the stationary Reynolds averaged 
Navier-Stokes (RANS) equations. We will show that the reduction of the system of equations in (9) or 
(8) in the infinite channels and pipes are the same (up to modified pressure) as the system of equations 
obtained in the case of NS-a (or the viscous Camassa-Holm equations (VCHE)), [6, 7, 8]. 
Let us begin by recalling the stationary RANS equations in channels and pipes (see, e.g., [38, 46]). We 
establish some notations: for a given function 4>(x,t) we denote by 

(0) (x) = fa) = lim 1 f 4>{x,t)dt (22) 

assuming that such a limit exists (see, e.g., [18] for the generalization of the notion of limit to make 
sense of infinite time averages.) The long (infinite) time average of the NSE, i.e. the stationary RANS 
equations, are given by 

(u ■ V)u = i/Au - Vp - (u - u) • V(u - u) 

V-ii = 0. {23) 

This averaging process yields the well known closure problem. The system above is not closed since we 
cannot express it solely in terms of u alone. The main idea behind turbulence modeling is to produce an 
approximate closed form for (23) in terms of u alone. 

3.1. The RANS equations for Turbulent Channel Flows. As might be expected from the visual 
appearance of the flow in experimental observations of turbulent Poiseuillc flows in infinite channel (sec, 
e.g., [38, 46]), the mean velocity in (23) for turbulent channel flows takes the form u = \U(z),0, 0] , 
where U(z) = U(—z), with mean pressure p = P(x,y,z). Using this classical observation, the RANS 
system (23) under such symmetry reduces to: 

-vU" + d z (wu) = -d x P 

8 Z (wv) = -d y P (24) 

d z (w 2 ) = -d z P 

where the prime (') denotes the derivative in the z-direction, and (u, v, w) T — u u is the fluctuation of 
the velocity in the infinite channel {(x, y, z) <G R, — d < z < d}. It is also observed from the experiments 
(see, e.g. [38, 46]), that the Reynolds stresses {wu) , (wv) and (w 2 ) are functions of the variable z alone. 
At the boundary, it is natural to impose the conditions U(±d) = (no-slip) and vU'(±d) = T T o, where 
To is the boundary shear stress. Using the boundary conditions (wu) (id) = (wv) (±<i) = 0, the Reynolds 
equations imply that (wv) = and P = Pq — r^x/d — (w 2 ) (z), with integration constant Pq. 
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3.2. The Reduced Simplified Bardinal Model for Channel Flows. For any turbulence model if 
an explicit analytical solution is available, then one can match this solution with the available physical 
experimental data to test its validity. Here we will show that the reduced simplified Bardina model under 
the channel symmetry admits the same exact equation as the reduced NS-a model. This is enough to 
show that the numerical solution of the reduced simplified Bardina model in the channel will match the 
experimental data for wide range of Reynolds number. For the simplified Bardina system of equations, 
under the channel symmetry, we denote by U the velocity u in (9) and we seek its steady state solutions 
in the form U = [C/(z),0, 0] , with even reflection symmetry condition U{z) = U{—z), and boundary 
condition U(±d) = 0. Under these conditions, the steady simplified Bardina equations reduces to: 

-vV" = -vU" + va 2 U"" = -d xP 

= -d y p (25) 
= -d zP 

where V = U — a 2 U" and p is a pressure function. Notice here that we need additional boundary condi- 
tions to determine V. Such boundary conditions are not yet available based on physical considerations. 
However, in this case, and under the symmetry of the channel, the missing boundary conditions come as 
free parameters that will be determined through a tuning process with empirical data. 

3.3. Identifying the Simplified Bardina Model with RANS - The Channel case. Following the 
idea of [6, 7, 8] we identify the systems (24) and (25) with each other, which is the essence of our closure 
assumption. We compare (24) and (25), and as a result, we identify the various counterparts as 

U = U 

d z {wu)=va 2 U"" + Pl 

d z {wv) = (26) 
V(P + (w 2 )) = V(p-pix) 
for some constant p\ . This identification gives 

{wv) = 0, 



2ttIII (^) 

— {wu) (z) = —p\z — va U 

and leaves (w 2 ) undetermined up to an arbitrary function of z. The identification in (26) is exactly 
the same (up to modified pressure and possibly (w 2 )) identification that was derived when identifying 
the NS-a model (VCHE) with the RANS equations in the channel symmetry in [6, 7, 8]. The same 
identification holds true in the case of the Leray-a model [10], the Clark-a model in [11] and the ML-a 
model in [26]. Therefore, similar to the earlier alpha models, the general solution of simplified Bardina 
and NS-a will be identical (up to a modified pressure) and in particular, the mean flows in both cases 
are the same functions. A similar result applies to turbulent pipe flows following the same argument and 
we will not include it here. For further details regarding the identification of the equations under the 
pipe symmetry, see [6, 7, 8, 26]. 

4. Existence and Uniqueness 

In this section we will prove the global existence and continuous dependence on initial data, (in partic- 
ular, the uniqueness of weak solution) of the system in (19). We will establish the estimates first for the 
finite dimensional Galcrkin approximation scheme and then using the appropriate Aubin compactness 
theorems (see for, e.g., [13, 44, 45]) we can pass to the limit. In this section, we fix T > to be arbitrarily 
large. 
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The finite dimensional Galerkin approximation, based on the eigenfunctions of the operator A, to (19) 
is : 

-^-(u m + a 2 Au m ) + vA(u m + a 2 Au m ) + P m B(u m ,u m ) = P m f 

u m (Q) = P m v? n . 

4.1. H 1 estimates. We take the inner product of the Galerkin approximation (28) with u m and use 
(17) to obtain 

\j t (\u m ? + a 2 |K„f) + v{\\u m f + a 2 \Au rn \ 2 ) = (P m f,u m ) = {f,P m u m ) = {f,u m ). (29) 
Notice that by the Cauchy-Schwarz inequality, we have 

lA-'fWAur, 



\{f,u m )\<[ | A -i/ 2/ ||| u m m || ( 3 °) 
and by Young's inequality we have 

\(f,U m )\<{ (31) 

We let K 1 = min { , — > , from the above inequalities we get 

v va z J 

j t {\u m \ 2 +a 2 \\u.J\ 2 )+iy(\\u m \\ 2 + a 2 \Au m \ 2 ) <K^ (32) 



^(K„| 2 + a 2 \\u m \\ 2 ) + u\ 1 (\u m \ 2 + a 2 \\u m \\ 2 ) < K\. (33) 



Applying Poincare inequality we get 

We then apply Gronwall's inequality to obtain 

\u m {t)\ 2 + a 2 \\u m (t)\\ 2 < e-^Qumm 2 + a 2 ||u m (0)|| 2 ) + ^(1 - e~" Alt ) (34) 

That is, 

\u m [t)\ 2 + a 2 |K(t)|| 2 < h := \u m \ 2 + a 2 \\u m \\ 2 + ^- (35) 

Thus, for t G [0, T], where T > arbitrary but finite, we get u m G L°°([0,T},V), where the bound is 
uniform in to, provided u m G V. 

4.2. H 2 estimates. Integrating (32) over the interval (t,t + r) for r > 0, we obtain 

nt+r 

(\\u m (s)\\ 2 + a 2 \Au m (s)\ 2 )ds < rK x + \u m (t)\ 2 + a 2 \\u m (t)\\ 2 



(36) 



< rK x + h- 

Now, take the inner product of the Galerkin approximation (28) with Au m to obtain 



~(\\u m \\ 2 + a 2 \Au m \ 2 ) + is{\Au m \ 2 + a 2 \A z ' 2 u m \ 2 ) + {B(u m ,u m ),Au m ) = (f,Au m ). (37) 



Notice that 
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Again by Young's inequality we have 

'^ 1/2/|2 ■ Via 3 / 2 * 



va 2 4 

\{f,Au m )\<{ (39) 



We denote by K 2 = min ^ , )■. Then we have 



~(\\u r J 2 + a 2 \Au m \ 2 ) + ^(\Au m \ 2 + a 2 \A^ 2 u rn \ 2 ) < K 2 + \(B(u m ,u m ), Au m )\. (40) 

Using Holder inequality, (11) and Young's inequality 

\(B(u ),Au m )\ < c\\ || 1/2 |Au m | 1/2 |Au m | 

= c||u m || 3 / 2 |A Wm | 3 / 2 (41) 
< ^\Au m \ 2 + c\\u m \\ 6 . 

Using the above estimates and (40) we obtain 

^(IKf + a 2 \Au m \ 2 ) + v(\Au m \ 2 + a 2 \A^ 2 u m \ 2 ) < 2K 2 + c\\u m \\ e . (42) 

We integrate the above equation over the interval {s,t) and use (35) and (36) to obtain: 

\\u m (t)\\ 2 + a 2 \Au m (t)\ 2 < \\u m {s)\\ 2 + a 2 \Au m (s)\ 2 + 2{t - s)K 2 + c(^j (t - s) (43) 
Now, we integrate with respect to s over the interval (0, t) and use (36) 

t{\\u m {t)\\ 2 + a 2 \Au m {t)\ 2 ) < I(fjjf 1 + k)+ t 2jjr 2 + ( /k > ) t- (44) 

for all t > 0. 

For t > we integrate (43) with respect to s over the interval (t — -^,t) 
-^-(\\u m (t)\\ 2 + a 2 \Au m (t)\ 2 ) 

< - (-k-Kx +k 1 S ] + 2K 2 (-^—] +c/ 



(45) 



v \v\\ ) \2v\\ ) \ a / \2^Ai 

Thus, from (44) and (45) we conclude: 

\\u m {t)\\ 2 + a 2 \Au m {t)\ 2 <k 2 (t) (46) 

for all t > 0. We note that, k 2 (t) enjoys the following properties: 

(i) k 2 (t) is finite for all t > 0. 

(ii) If u m e V, but u in i D(A). then the lim t ^ + k 2 (t) = oo. 
(hi) limsup^oo k 2 (t) < oo. 

Remark 1. ^From (43), one can observe that if u m S D(A), then u m (-) is bounded uniformly in the 
L°°([0,T];D(A)) norm, independently of m. On the other hand, if u m e V, but u m <£ D(A), we 
conclude from the above that u m £ L% c ((0, T],D(A)) n L 2 ([0, T], £>(A)). 
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In order to extract convergent subsequence by using Aubin's lemma (see [13, 34, 44]), we need to establish 

dv m du m 
estimates tor — ; — , — ; — . 

dt dt 

dv 

—j^- = -vAv m ~ B(u m , u m ) + P m f (47) 

Take the D(A)' (the dual of the space D(A)) action of the equation above with it) 6 D(A), we observe 
that 

|(P TO j»| = \(f,P m w)\ < lA-'fWAw] < \^\f\\Aw\ = L 2 \f\\Aw\ (48) 
and using (14), we have 

\(P m B(u m ,u m ),w)\ < c\u m \ 1/2 \\u m \\ 1/2 \\u m \\\\w\\ 



c u m ' u 



Vaii,. i|3/2| 



w\\ (49) 



= C A- 1/2 | Mm |V2|| Mm ||3/2|^|. 

By (35), ||um||L°°([o,T];V) is bounded uniformly with respect to m. Thus by (49), wc can deduce that 
\\P m B(u m ,u m )\\L 2 ([o,T];D(Ay) is a l so bounded uniformly with respect to m. Now, the uniform in m, 
L 2 ([0,T]; D(A)) bound for u m implies that ||^m||L 2 ([o,T];ff) is uniformly bounded, which in turn implies 

that ||^4wm||L 2 ([o,T]:D(A)') is uniformly bounded, as well. Thus, we conclude, || ^ \\l 2 {[o,t];D{A)')i an d i n 

particular, || -j- \\l 2 ([v.t];H)- ar c uniformly bounded with respect to m. By Aubin compactness theorem 
(see, e.g., [13, 34, 44]) we conclude that there is a subsequence u m '(t) and a function u{t) such that 

u m >(t) -> u(t) weakly in L 2 ([Q, T];D(A)) 

-» «(*) strongly in L 2 ([0, T];V) (50) 

tw-uin C([0,T];H), 

or equivalently, 

v m ,{t) -> v(t) weakly in L 2 ([0, T];H) 

v m ,(t) -> v(t) strongly in L 2 ([0, T]; V) (51) 
iW -^vin C([0,T];D(A)'). 
We relabel u m < and v m i with u m and v m respectively. Let w £ -D(^4), then we have 



(v m (t),w)+v {v m {s),Aw)ds + / (i?(u m (s),u m (s)),P m ti;)ds = {v m (t ),w) + (f,P m w)(t-t ) (52) 
for all t ,t E [0, T]. The sequence u m (i) converges weakly in £ 2 ([0, T]; if) and thus, 

lim f (v m (s),Aw)ds = [ (v(s),Aw)ds. (53) 



m — >oo 



Also, by (51), v m (t) converging weakly in L 2 ([0,T]; H) implies that there is a subsequence of v m , which 
we relabel as v m , which converges a.e. t G [0, T] to v(t) in H' ~ ff. Thus, we conclude that 

(«,„(*), ttf) (u(i), to), and 

(v m (t ),w) -» (u(io),u>) 

for a.e. t,to G [0, T]. On the other hand, 

t 

(B(u m (s),u m (s)),P m w) - (B(u(s),u( S )),w)d S \ < /&> + J« +4 3 ). (55) 
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Using (18) and Agmon's inequality (12), we get 



(B(u m (s),u m (s)),P m w - w) ds 



t 



< c 



< c 



(s)|||u m (s)||||P m w - w\\ L oa(Q)dS 



(56) 



ft \ 1/2 / ,t \ 1/2 

f |w m (s)| 2 ds / \\u m (s)\\ 2 ds) \P m w~w\ 1 ' i \A{P m w - w)\ 3/i 



Since u m is bounded uniformly in L°°([0,T];V) independent of m , and thus is bounded uniformly in 
L°°([0, T];H) thanks to Poincare inequality, we get linim^oo 1$ = 0. 

Again, using (18), Agmon's inequality (12), and Poincare inequality, we get 



/(2) 



(B(u m (s) ~ u(s),u m (s)),w) ds 



to 



< 



c / \u m {s) - u m (s)\\\u m (s)\\\\w\\ L ^^ n) ds 



to 



(57) 



< c 



\u m (s) - u(s)\ 2 ds 



1/2 



to 



1/2 



\\u m ( S )\\ 2 ds) \-V*\Aw\. 



Now since u m — > u strongly in L 2 ([0,T];V) (thus in L 2 ([Q,T]; H) ) and u m is bounded uniformly 
independent of m in L°°([0, T]; V) , we get that lim„ 



j(a) _ a 



j(3) 



(B(u(s), u m (s) — m(s)), w) (is 



to 



to 



<c / |u(s)|||u m (s) - u{s)\\\\w\\ L ^(n)ds 

1/2 / 

|u(s)| 2 rfs ] 

'to 



< C 



\ 1/2 ft 1 \ 1/2 

U ( S )| 2 rf S J (jf \\u m (s)-u(s)\\ 2 ds) \; 1/A \Aw\ 



where we applied (14) and (18) in the second and third inequality, respectively. Now, u Tl 

in L 2 {[0, T];V) implies that liim^^ 1$ = 0. 

^From the above calculations, we have that for a.e. to, t £ [0, T], 



(58) 



u strongly 



(v(t), w) — (v(to), w) + v I (v,Aw)ds+ / (B(u(s), u(s)), w) ds 



(f,w)ds. 



(59) 



for every w £ D(A). To show that v £ C([0,T]; V), and hence, u £ C([0,T];V), we want to show that 



the viscous term v J (v(s), Aw)ds and the nonlinear term J (B(u,u) 



v I (v(s),Aw)ds 



< V 



t \ V2 

2 



Ks)| 2 ds 



to 

t \ 1/2 

2 



, ds — ► as t — * to. 



\Aw\ 2 ds 



ID(A) 



— > as t —>to, 



(60) 



since w € L 2 ([0,T];iT) and w G 



((B(w( s ), u (s)), W )ds) <Hi-(n) / |«(a)rds 



1/2 



t \ 1/2 

2, 



|u(s)|| 2 <is 



as t^t , (61) 



since u G L°°([0,T];V). Thus, this implies that for a.e. t £ [0,T], (v(t),w) — > (u(i ), as t — > i , for 
every w e In particular, £ H C V and to e -D(A) C V, implies that for a.e. t £ [0,T], 
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(v(t),w) v , — ► (v(to),w) v , as £ — > to, for every w G D{A). Since -D(-A) is dense in V, for any test 
function (j> £ V and for every e > 0, there exists a u> 6 -D(A) such that — (f>\\ < e/(M + 1), where 
M = 2sup te r 0)T i ||u(£)||y. Thus for every 4> G V 

(v(t) - v(to),(f>) v , | < | (v{t) - v(t ),w) v , | + | (v(t) - v(t ), w - <j))y, \. (62) 

The first term goes to zero as t — > to since w G -D(^l). For the second term, we have 

| (v(t) - v{t Q ),w - <j>) v , | < ||w(*)-t;(to)||v/||«;-0|| < M\\w - 0|| < c. (63) 

Since e > is arbitrary, we conclude that (f(t) — ^(to) ; y — * 0, as t — > to: for a ll S V. Hence, 
« £ C([0,T]; V) and in particular, u £ C([0,T];F). 

To summarize: we have established above the global existence of weak solution of the simplified Bardina 
system by the standard Galerkin approximation scheme together with some useful a priori estimates. 

Theorem 3. (Global existence and uniqueness) Let f S H and u m G V. Then for any T > 0, (19) 
has a unique weak solution u in [0,T]. 

To complete the proof of the theorem above, we are left to prove the uniqueness of weak solutions. 

Uniqueness of Weak Solution. Next we will show the continuous dependence of the weak solutions 
in the appropriate norm specified below, on the initial data, and in particular, the uniqueness of weak 
solutions. 

Let u and u be any two weak solutions of (19) on the interval [0, T], with initial values u(0) = u ln G V 
and u(0) = u m G V, respectively. Let us denote by v = (u + a 2 Au), v — (u + a 2 Au), Su = u — u, and 
by Sv — v — v. Then from (19) we get: 

j t 5v + uASv + B(Su, u) + B(u, Su) = (64) 
By taking the D(A)' action of (64) with Su, 

/■^-5v 7 5u\ + v (ASv,5u) D{A y + (B(Su,u),5u) + (B(u,Su),Su) = (65) 

\™ / D(A)' 

and by applying a Lemma of Lions-Magenes concerning the derivative of functions with values in Banach 
space, (cf. Chap. III-p.l69-[44]) and by the property of the bilinear form B, (17), we get: 

\ji\H 2 + a 2 \\Suf) + KIIHI 2 + a 2 \ASu\ 2 ) + (B(6u, u), Su) = 0. (66) 
Dropping the nonnegative viscous term, and by using property (14), we get 



i|(IH 2 + « 2 IIHI 2 )< C |IHI 2 H 1/2 IHr /2 



< cA^ 1/4 ||5u|| 2 | 



(67) 



By Gronwall inequality, we obtain: 



(\Su(t)\ 2 + a 2 \\Su(t)\\ 2 )<(\5u(0)\ 2 + a 2 \\Su(0)\\ 2 )exp[ I C ^ u( f^ ds) . (68) 



In (68), since u G L°°([0, T];V), we have shown the continuous dependence of the weak solutions on the 
initial data in the L°°([0,T]; V) norm. In particular, in the case we have the same initial data, we have 
||<5w(i)|| 2 = 0, which implies that we have u(t) = u(t), for all t G [0, T]. 
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5. Global Attractors, Their Dimensions and Connection to Dissipation Length Scales 

Now that we have established the global wcll-poscdness to the simplified Bardina model, in this section 
we will show the existence of global attractor A C V for the system (19), its finite Hausdorff and fractal 
dimensions, and the physical relevance of this finite dimension of global attractor to the concept of "finite 
dimensionality" of turbulent flows. 

Following standard techniques, the method that we will use to estimate the dimension of the global 
attractor stems from the following lemmas (see [13, 32, 45] and [16], respectively): 

Lemma 4. (The Lieb-Thirring inequality). Let {"0j j"j=i ^ e an orthonormal set of functions in 
(H) k — H © H ■ ■ ■ © H . Then there exists a constant Clt , which depends on k, but is independent of 

k-times 
N , such that 

)5/3 
N 
dx < G LT V / (Vi/>j(x) : Vipj(x))dx. (69) 

Lemma 5. Let {<^j}|Li & V be an orthonormal set of functions with respect to the inner product [•, •]: 

Let ipj(x) = ((f) j (x) , a , a , a ) , and 4> 2 (x) = X)jLi(^j(^) • <t>j(p))- Then there exists a 
constant Cf, which is independent of N , such that 

(N \ 1//2 

^Jjy^{x):V^{x))dx\ . (70) 

To start the study of the finite dimensionality of the Hausdorff and fractal dimensions of the global 
attractor, first we recall that, from the existence and uniqueness properties of the solutions to (19), we 
get a semi-group of solution operators, denoted as {S(t)}t>o, which associates, to each u™ s V, the 
semi-flow for time t > : S(t)u m = u(t). We are now ready state and prove the following theorem: 

Theorem 6. There is a compact global attractor A C V, in terms of the solution u, for the system (19). 
Moreover, we have an upper bound for the Hausdorff and fractal dimension of the attractor A 

d H (A) < dp (A) < c G 6 / 5 (—1—) = c G 6 / 5 (-) ' (71) 



a 



l/l 

where G = tttt is the Grashoff number. 

2\ 3 / 4 

Proof. The first requirement to show the existence of the nonempty compact attractor is to show that 
we have an absorbing ball in V and D(A) and that the semigroup {S(t)} t >o defined above is compact 
(see, e.g., [13, 22, 39, 42, 45]). This can be established from the previous a priori estimates. First, let 
us show that there is an absorbing ball in V and D(A). By (34), (50), and the fact that |u m (0)| < |u(0)| 
and ||u m (0)|| < ||m(0)|| we have, by passing to the limit with m — > oo, 

KOI 2 + a 2 \\u(t)\\ 2 < e-" Alt (K0)| 2 + a 2 |K0)|| 2 ) + ^(1 ~ e"^ 1 *). (72) 



Choose t large enough such that e 1/Ai *(|m(0)| 2 + a 2 ||u(0)|| 2 ) < — then we have 

vXi 

\u(t)\ 2 + a 2 \\u(t)\\ 2 <2^, (73) 



14 



Y. CAO, E. LUNASIN, AND E.S. TITI 



\A-ifi \A-vy\* \ 

where wc recall K\ — mm < 7. — , > ■ In particular. 



limsup (\u(t)\ 2 + a 2 \\u{t)\\ 2 ) < 2^- =: R 2 V . (74) 

Therefore, the system (9) has the ball By(0) in V of radius Ry as an absorbing ball in V. 

Proving the existence of absorbing ball Bd(a)(0) m D{A) is similar. By (45) and (46) we conclude that 

limsup {\\u(t)\\ 2 + a 2 \Au{t)\ 2 ) < limsup k 2 {t) =: R 2 D{A) < 00, (75) 

t — >oo t — >oo 

and therefore we have the ball Bd(A)(0) hi D(A) with radius Rd(A) as an absorbing ball in D(A). 

Now applying Rcllich lemma [1] wc have that S(t) : V — > T>(A) CC V, for t > 0, is a compact 

semigroup from V to itself. What is left is to show that indeed we have a nonempty compact attractor. 

v 

Since S(t)Bv(0) C Bv(0), it follows that for each s > the set C s := Dt> s S(t)Bv{0) is nonempty and 
compact in V. By monotonicity of C s for s > 0, and by the finite intersection property of compact sets, 
we see that 

^=HU S(t)B v (0) V c V (76) 

s>0 t>s 

is a nonempty compact set in V and indeed is the unique global attractor in V. 

We are now ready to give an upper bound estimate to the Hausdorff and fractal dimensions of the 
global attractor. As mentioned above, we will use the trace formula (sec, e.g., [12, 13, 45]) to establish 
this estimate. 

The first step in this estimation is to do linearization about a solution. We note that in order to 
apply the techniques in [12, 13, 16, 45], we need that the mapping S(t) : V — > V is differentiable with 
respect to initial data. Following similar ideas of energy estimates in the proof of uniqueness of weak 
solutions in the previous section, one can show that S(t)u m is differentiable with respect to u m , when 
u m e A. Thus said, we linearize the viscous simplified Bardina model (19) about a solution u{t) (or 
v(t) = u(t) + a 2 Au(t)) 

^-5v + uA5v + B(Su, u) + B(u 7 Su) = 

dt (77) 
Sv(0) = Sv m = 5u ln + a 2 A8u m 

where 6v is a perturbation satisfying (77) and is given by Sv = 8u + a 2 ASu. With this relationship, Su 
evolves according to the equation 

^-Su + vASu + (I + a 2 A)- 1 [B{5u, u) + B(u, Su)} = 

at \J 8) 



Su(0) = 5u l 



which we write symbolically as 



4<5w + T(t)Su = , 

dt v ; (79) 

<Ju(0) = Su m , 

where T(t)ip = vAip + (I + a 2 A)- 1 [B(ijj, u(t)) + B{u{t),ip)]. Let Sui(0), j = 1, . . . , N be a set of 

linearly independent vectors in V and let Suj(t) be the corresponding solutions of (78) with initial value 
8uj(0) for j = l,...,N. Let 

T N (t) = Tracc(P w (i) o T(t) o Pjv(i)) (80) 
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where P/v(i) is the orthogonal projection of V onto the span {6vi(t), Sv2(t), . . . , SvN(t)}. We shall denote 

by {<fij(t)}j=i JVj an orthonormal basis, with respect to inner product [•,•] = (•,•) + ct 2 ((-,-)) of the 

space PnV = span{Svi(t), Sv 2 (t)}. ^From (80) we have 



AT 



3=1 

iV 

= £ v{A<f>j , (j)^] + [(I + crA^B^^u),^ + [(J + a 2 ^)- 1 ^, h 



3 1 (81) 
= "Yl IM'.^'] + E ( fl fe> u )^i) + E (fl(u,&),&) 

3=1 3=1 3=1 

AT AT 

3=1 3=1 

By the definition of the inner product [■, ■], we have 

N N N N 

[MM = E (MM +c?Y^(Mi,Mj) = E / : ViJj(x,t))dx =: Q N (t) (82) 

3=1 3=1 3=1 3=1 Jn 

where, 

( , d , d d 

i'j = I fa, «; 

Note also that 
Setting 

we have 



*> = V ** a ~^ h a d^ h a ^ ' • (83) 



(i>i,'M=S j k. (84) 



AT 
3 = 1 



3^(t) = i/Qj>r(t)+Ww(t). (85) 
We denote by ip 2 := YljLi V*j " V'j- F° r 7£/v(*) we have 

AT » JV „ AT 



3=1 Jn 3=1 ^ 3=1 



(86) 



By the estimates so obtained above we finally find 



cum 



TArW>^Qiv(*)-^|Kt)|| 2 . (87) 
Z Ivor 

By the asymptotic behavior of the eigenvalues of the operator A (see (10)) and (84) we get 

AT N 

Qn® = E II 2 > E A ^ coAiTV 5 / 3 . (88) 
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Now, by the trace formula (see, e.g., [12, 13, 45] and the references therein) if N is large enough so that 

1 f T 

liminf- / T N {t) dt > (89) 

then N is an upper bound for the Hausdorff and fractal dimensions [13, 45], (see also [5]), of the global 
attractor. 

Thus, by (87) and (88) it is sufficient to require N to be large enough such that 

^AxiV 5/3 > sup limsupi / ^f^\\u{t)\\ 2 dt. (90) 



On the other hand, using Holder inequality we get from (36) 

^rs^'^i^^ (9i » 

which implies 

¥ 5 /3 > c\m\ 2 i/i 2 cm ^r 2 

- l/*A?0« ^4 A 3/2 A 3/2 a6 - A 3/2 a6 

and we recall that |fi| = (2irL) 3 and that Ai = L~ 2 , thus, 

( 

A~fV 8 /5 



(92) 



N > G 6/5 -^— (93) 



;,From this wc deduce that 

A„( A\ <r A„( A\ <r n®/'_ 

A? /5 a 18 / 5 



d H (A) < d F (A) < G 6/5 , 9A , C o;c ■ (94) 

□ 



The interpretation of the upper bound estimate that we get for the Hausdorff and fractal dimension of 
the global attractor in terms of small scales is important in showing the finite dimensionality of flows and 
in particular in showing the numerical computability of the turbulence model. To do this, wc interpret 
the estimate for the attractor dimension in terms of the mean rate of energy dissipation of the simplified 
Bardina model. Following [16] wc define the corresponding mean rate of dissipation of "energy" for the 
simplified Bardina model (see (29)) as 

i r T 

e = L~ 3 v sup limsup- / (\\u(s)\\ 2 + a 2 \Au(s)\ 2 )ds. (95) 



T 



Thus, and in analogy with the Kolmogorov dissipation length in the classical theory of turbulence, we 
set the dissipation length scale for the simplified Bardina model as 

h = (j) ■ (96) 

Identifying the dimension of global attractor with the number of degrees of freedom, we will show that 
the number of degrees of freedom for the simplified Bardina model is bounded from above by a quantity 
which scales like (L/a) 12 / 5 {L/l d ) 12 / 5 . 

In fact, in view of (95) we can write (91) as follows 



(97) 
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Using this in (90) and recalling (96) we obtain the following estimate for the dimension of the global 
attractor and hence the upper bound on the number of degrees of freedom in the simplified Bardina 
model is: 

/r\ 12/5 12/5 

d H (A)<d F (A)<c(-) f-J . (98) 

We remark that following the recent work of [21] one can also interpret this bound in terms of the 
Reynolds number. 



6. Energy Spectra 

Turbulent flows are characterized by the presence of wide range of eddy sizes starting from the size 
of the flow domain, say 2irL in our case, to much smaller scales, which become progressively smaller 
relative to 2itL as we increase the Reynolds number. It is important to examine how the energy in 
a turbulent flow is distributed among these different size eddies by considering the energy spectrum. 
Following similar arguments to those presented in [15] and [17] (see also [10, 11, 18, 26]) we will study 
in this section the energy spectra of the simplified Bardina model. We will obtain our results about the 
decay of the energy spectrum for the filtered velocity u following similar techniques as those in the NS-a 
[17]. In particular, we observe that there are two different power laws for the energy cascade. For wave 
numbers k <C 1/a, we obtain the usual fc~ 5 / 3 Kolmogorov power law. This implies that the large scale 
statistics of the flow, in particular for those eddies of size greater than the length scale a, are computed 
consistent with the Kolmogorov theory for 3D turbulent flows. On the other hand, for k ^ 1/a, that is 
for eddies smaller than the length scale a, we obtain a steeper power law. The steeper spectral slope for 
wave numbers k 3> 1/a implies a faster decay of energy in comparison to DNS, which suggests, in terms 
of numerical simulation, a smaller resolution requirement in computing turbulent flows. For this reason, 
we suggest that the simplified Bardina model is a good candidate for a subgrid scale model of large eddy 
simulation of turbulence. To start, we define some notations: 

6(u,w,?ii) = (£?(«, i>), u>), 

' ' - -ik-. 



(2-kLY 



Q 



(2ttL)2 



k<\j\<2k 



Vie 



k<\j\<2k 

U k =E U J' V k =J2 v j 

j<k j<k 

u k = Yl u i> v k = E v i- 

2k<j 2k<j 

There are three flow regimes that we need to consider to analyze the energy spectra. These are the 
flow regimes where energy is produced, where energy cascades (i.e. inertial range) and, where energy 
dissipates and decays exponentially fast (i.e. dissipation range). We split the flow into three parts 
according to the three length scale ranges. Assume kf < k, where kf is the largest wavenumber involved 
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in the forcing term. Thus, 



U = uf + U k + Uu 



v = v< + v k +v>. 

The energy balance equation for the simplified Bardina model for an eddy of size k~ x is given by 

\ ^ Vk,Uk ^ + ^(~ Aufe ' Ufc ) = T k~ T 2k, (99) 

where, 

T k := -b(u^,u^,u k ) + b(u k + u^u k + u k ,u k ). (100) 



We can interpret T k as representing the net amount of energy per unit time that is transferred into 
wavenumbers larger than or equal to k. Similarly, T2 k represents the net amount of energy per unit 
time that is transferred into wavenumbers larger than or equal to 2k. From these definitions T k — T 2k 
represents the net amount of energy per unit time that is transferred into wavenumbers between [k, 2k). 
Taking an ensemble average (long time average) of (99) we get: 

v ({-Av k , u k )) = (T fc ) - (T 2k ) . (101) 

We define the energy of eddy of size 1/k as 

E a (k) = (l + a 2 \k\ 2 ) J2 I"/- 
Ul=fc 

This definition arose from the fact that we consider \u\ 2 + a 2 ||w|| 2 as the "energy", since this is the 
conserved quantity in the simplified Bardina model equation (see section 7). Using this definition, we 
can now rewrite the time-averaged energy transfer equation (101) as 

2k 

2 



uk A E a {k) ~ v J k z E a (k)dk ~ (T fe ) - (T 2k ) . 

Thus as long as vk 3 E a (k) << (T k ) (that is, (T^fe) ~ (Tk), there is no leakage of energy due to dissipation), 
the wavenumber k belongs to the inertial range. Similar to the other alpha subgrid scale models, it is 
not known what is the correct averaged velocity of an eddy of length size That is, we do not know 
a priori in these models the exact eddy turn over time of an eddy of size k^ 1 . As we will see below, we 
have a few candidates for such an averaged velocity. Namely, 

\ 1/2 

^ = (73 / l^l 2 ^) ~(/ {l + o?k 2 )E a {k)\ ~(k(l + a 2 k 2 )E a (k)) 1/2 , 

1/2 

U l = (-k I u *- v * dx ) ~ ( / E «( k ) ) - (kE a (k)) 1/2 , 




kE a {k) 
- a 2 k 2 



1/2 



that is, 



^- (1 + ^^ (n = M ' 2) - (102) 
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In the inertial range, the Kraichnan energy cascade mechanism states that the corresponding turn over 
time of eddies of spatial size 1/k with given average velocity as above is about 



1 



(l+a 2 ^ 2 )^- 1 )/ 2 



kur 



(n = 0,l,2). 



Therefore the energy dissipation rate e is 



1 f 2k 

— / E a (k)dk 

T k Jk 



k 3 / 2 (E a {k)y/ 2 

k^ 2 (E a (k)f 2 



(1 



and hence 



E a {k) 



a 2 k 2 )^ 



-l)/2' 



(103) 



^(l + a 2 /^"- 1 )/ 
W 3 



Note that the kinetic energy spectrum of the variable u is given by 

2/3 

E a (k) 



E u (k) 



£5/3 ^ 



1 + a 2 ^ 



2/3 



when fa « 1 , 
when fca>l. 



Q ,2(4-ra)/3£(13-2n)/3 ' 

Therefore, depending on the appropriate average velocity on an eddy of size for the simplified 
Bardina model, we would get the corresponding energy spectra which has a much faster decaying power 
law fe( 2n - 13 )/ 3 , (n = 0, 1, 2) than the usual Kolmo gorov k 5 / 3 power law, in the subrange to > 1. This 
signifies that the simplified Bardina model, like the other alpha models, is a good candidate subgrid scale 
model of turbulence. 



7. Global Existence and Uniqueness of the Inviscid Simplified Bardina Model 

In this section, we will established the global existence and uniqueness of the inviscid simplified Bardina 
model using the classical Picard iteration method. The inviscid simplified Bardina modle is equivalent 
to the functional differential equation 

dv 

— + B(u,u) = f, 

v = u + a 2 Au, ( 104 ) 
v{0) = v in =u ln + a 2 Au m , 
where, for simplicity, we assumed / to be time independent. 

Theorem 7. (Short time existence and uniqueness). Let v m G V', and f £ V . There exists a 
short time T*(\\v m \\v') such that the equation (104) has a unique solution ufC 1 ([— T», T*], V), that is, 

Proof. We will use the classical Picard iteration principle (see, e.g., [41]) to prove the short time existence 
and uniqueness theorem. Namely, it is enough to show that the vector field N(v) = f — B(u, u) is locally 
Lipschitz in the Hilbert space V . From the classical theory of ordinary differential equations we consider 
the equivalent equation for (104) 

v(t)=v in - f B{u{s),u{s))ds + f t. (105) 
Jo 

Notice that v € V implies that u £ V and thus by Lemma 1 B(u, u) € V. As a result the equation 
above makes sense in the space V . Let vi,v 2 £ V, and consequently ui,u 2 £ V. By (14) and Poincare 
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IliVM - N(v a )\\ v > < -nxll u i ~ "all < ^Iki - v 2 \\ v ,. (107) 



inequality, we have 

\\N(vx)-N(v2)\\v> = \\B( Ul , Ul ) - B(u 2 ,u 2 )\\v< 

= SUp I (B(ui - U2,U 2 ) + B(U!,UI - U2),U)) V , I 

{wev,\\w\\=i} (106) 

< ^4ll«l-«2||(ll«l|| + ll«2||)- 

For any large enough R such that [|ui||, ||w2|| < R, we have 

4ci? n „ ci?, 

Here we used the fact that \\v\\v is equivalent to ||u||. Equation (107) implies that N(v) is locally Lips- 
chitz continuous function in the Hilbert space V . By the classical theory of ordinary differential equa- 
tions, the equation (105) has a unique fixed point in a small interval [— T*,T*] and v G C([— T*, T m ], V) 
(see, e.g., [41]). In particular, since B(u(s),u(s)) is a continuous function with values in V and the forc- 
ing / assumed to be time independent, equation (105) implies that the left hand side v(t) is differentiable 
and 

nil 

^ = -B(u, u) + f 

dt y ' 3 (108) 

v{0) =v in . 

This implies the local-in-time existence and uniqueness of solution v G C 1 ([— T t , T«], V), and hence, 
u G C 1 ([— T*, T*], y), to the inviscid simplified Bardina model (104) or (108). We will next show that, 
in fact, we have global existence. To show global existence to (104) or (108) it is enough to show that 
on the maximal interval of existence, ||i>(i)l|y remains finite. Let [0,T max ) be the maximal interval of 
existence. If T max = +oo, then there is nothing to prove. Suppose, for the purpose of contradiction, that 

T max < oo. (109) 

This implies that limsup t ^ T - II^C*) II V = oo. By the equivalence of the norms ||v(t)||y and ||u(i)||, we 
conclude that also 

limsup ||u(i)|| = oo. (HO) 
We will derive a contradiction to the conclusion in (110). 

Notice that on [0, T max ), u G C([0, T max ), V), hence we can take the action of (104) or (108) on u(t). We 
get, by (17) 

(c^' U ) = ~ ( B ( u > u )> u )v> + (f> u )v = (f> u )v> ( m ) 

Thus, we have 

^(M 2 + a 2 IMI 2 )<||/||HM|. (H2) 

Let eo be a positive constant which has the same units as |m| 2 . From (112) we have 



i!(M 2 + a 2 || u || 2 + eo )<||/IMMI 



< 



(113) 
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Denote by z 2 := \u\ 2 + a 2 ||tt|| 2 + eo- Then we can rewrite (113) as 

dz 
~di 



dz < Jj/jjy^ (n4) 



Consequently, 



z(t) < z(0) + MMLt, (115) 



n 



for all t < T max . Therefore, by letting e — > we obtain 



and in particular, 



This implies that 



KOI 2 + a 2 \\u(t)\\ 2 < \u(0)\ 2 + a 2 \\u(i))\\ + MML t , (n 6 ) 



\u(t)\ 2 + a 2 \\u(t)\\ 2 < \u(0)\ 2 + a 2 \\u(0)\\ + ^^T max =: K. (117) 



limsupKOI + a 2 \\u(t)\\ 2 < K. (118) 

This is a contradiction to the conclusion (110). □ 
To summarize, we have established the proof to the following theorem: 

Theorem 8. (Global existence and uniqueness) Let f G V' and v m G V . Then the system in 
(104) has a unique solution v G C 1 ((— 00, 00), V) ( or equivalently, u G C 1 ((— 00, 00), V)). 

We observe that the inviscid Bardina model, (104), is equivalent to the following modification of the 
3D Euler equations 

-a 2 A— + — + (u ■ V)« + Vp = f 

V-u = (119) 
u(x,0) = u m . 

In particular, it is equal to the Euler equations when a — 0. Therefore, wc propose the inviscid simplified 
Bardina model as regularization of the 3D Euler equations that could be implemented in numerical 
computations of three dimensional inviscid flows. The analytical study of the regularity of the solutions 
of the inviscid simplified Bardina model, and in particular the limit of its solutions, as a — > 0, to the 
solutions of the Euler equations will be reported in a forthcoming paper. 

Inspired by the above model, (sec also [27, 37]), we propose the following regularization of the 3D 
Navier-Stokes equations 

-a 2 A— + — - i/Au + (u ■ V)u + Vp = f 
at at 

V-u = (120) 

u(x, 0) = u m 

subject to either periodic boundary condition or the no-slip Dirichlct boundary condition u\qq, = 0. In the 
presence of physical boundaries the above regularization (120) of the Navier-Stokes equations is different 
in nature from the hyperviscosity regularization of Lions [33] , or any of the other alpha regularization 
models, because it does not require any additional boundary conditions. It is also simpler than the 
nonlinear viscosity model of Ladyzhcnskaya [29, 30] and Smogarinsky [43]. We will study the analytical 
and long-term properties of (120) in a forthcoming paper. 
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